Calculation of the connective constant for self-avoiding 
walks via the pivot algorithm 

Nathan Clisby 

ARC Centre of Excellence for Mathematics and Statistics of Complex Systems, 
Department of Mathematics and Statistics, 
The University of Melbourne, VIC 3010, Australia 

February 8, 2013 



Abstract 

We calculate the connective constant for self-avoiding walks on the simple cubic lattice to unprece- 
dented accuracy, using a novel appUcation of the pivot algorithm. We estimate that = 4.684 039 931± 
0.000 000 027. Our method also provides accurate estimates of the number of self-avoiding walks, 

even for walks with millions of steps. 

Keywords self-avoiding walk; connective constant; Monte Carlo; pivot algorithm; approximate enu- 
meration 



1 Introduction 

The self-avoiding walk (SAW) on a regular lattice is an important model in statistical mechanics with 
a long history [1]. An A^-step SAW is a map uj from the integers {0, 1, • • • , A^} to sites on the lattice, 
with w(0) conventionally at the origin, + 1) — u:{i)\ — 1, and u}{i) ^ LLi(j) V« 7^ j. SAW is a 
topic of much current interest: see [2] for a recent review of rigorous results, and [3] for an overview 
of self-avoiding polygons (SAP) which has broader scope, including numerical aspects of SAP and to a 
lesser extent SAW. 

The most important quantities which characterize SAW are the number of SAW of length N, cm, and 
measures of the size of the walk, such as the square end-to-end distance. The asymptotic behavior of cjv 
on the simple cubic lattice is believed to be 

CNr^A^Ji''N^-^{l + 0{N-'^^)), (1) 

where the connective constant n and amplitude A are lattice dependent, the critical exponent 7 is uni- 
versal, and Ai is the exponent of the leading correction to scaling. There are also sub-leading analytic 
corrections to scaling, and a contribution from the so-called anti-ferromagnetic singularity; see for ex- 
ample [4] for more details on the asymptotic form of cjv- 

Enumeration is a particularly powerful method for studying SAW on two-dimensional lattices, where 
the finite lattice method is highly effective [5, 6, 7]. The best estimate for ji on the square lattice comes 
from enumerations of self-avoiding polygons to 130 steps [8], leading to the highly accurate estimate 
ji = 2.638 158 530 35(2). For the simple cubic lattice, the best estimate for /i comes from PERM Monte 
Carlo simulations [9]: /z = 4.684038 6(11). The most powerful known enumeration method for three- 
dimensional lattices is the length-doubling algorithm [10], which combines brute force enumeration with 
the inclusion-exclusion principle in a novel way. SAW on the simple cubic lattice have been enumerated 
to 36 steps, with C36 = 2 941 370 856 334 701 726 560 670, and /x = 4.684 040 1(50) [10]. 

In this paper we will obtain a highly accurate estimate of fj, for SAW on the simple cubic lattice using 
a Monte Carlo algorithm. Our method can also be used to estimate the number of self-avoiding walks. 
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2 Method 



Our method to calculate ji for SAW combines four key ideas: 

1. Use of the pivot algorithm, the most powerful known method for sampling SAW; 

2. A novel computer experiment which involves a telescoping sum that eliminates corrections to 



3. The adoption of scale-free moves to efficiently calculate the observable of interest; 

4. Partitioning CPU time between different sub-problems in an optimal way. 
We now describe each of these aspects in turn. 

2.1 The pivot algorithm 

The pivot algorithm is an extremely powerful method for sampling SAW in the canonical ensemble. 
It was invented by Lai [11], but the true power of the method was only appreciated after the ground- 
breaking work of Madras and Sokal [12]. Recently, the implementation of the pivot algorithm has been 
improved to make it even more powerful [13, 14, 15]. The recent improvements make it an extremely 
attractive prospect to utilize the pivot algorithm whenever possible. 

The pivot algorithm is a Markov chain Monte Carlo algorithm which works in the set of self-avoiding 
walks of fixed length, where the elementary move is a pivot as described below. The pivot algorithm 
generates a correlated sequence of SAW via the following process: 

1. Select a pivot site of the current SAW according to some prescription - usually uniformly at ran- 



2. Randomly choose a lattice symmetry (rotation or reflection); 

3. Apply this synmietry to one of the two sub- walks created by spUtting the walk at the pivot site; 

4. If the resulting walk is self-avoiding: accept the pivot and update the configuration; 

5. If the resulting walk is not self-avoiding: reject the pivot and keep the old configuration; 



The pivot algorithm is ergodic, and satisfies the detailed balance condition which ensures that SAW are 
sampled uniformly at random [12]. 

After a successful pivot, global observables, such as the square end-to-end distance, change signifi- 
cantly and are essentially uncorrelated. This observation is equivalent to the statement that the integrated 
autocorrelation time for a global observable A, Tint(^)> is of the same order as the mean time for a suc- 
cessful pivot. In the language of [12], once a successful pivot has been made the resulting configuration 
is "essentially new" with respect to global observables. For SAW on the simple cubic lattice the proba- 
bility of a pivot attempt being successful is 0{N~p), withp w 0.11. Therefore global observables have 
Tint = 0{N^); see [12] for extensive discussion. 

For local observables, such as the angle between the 37th and 38th steps of a walk, one may need 
0{N) successful pivots before the observable changes. Consequently Tint = 0{N^~^p) for local observ- 
ables. 

2.2 Telescoping observable 

Given wiilks uii and uj2, we define a concatenation operation by placing the root point of oji at the origin, 
and the root point of at (1, 0, 0). We denote the resulting walk as uji o uj2- Under this definition of 
concatenation, walks of M and A'^ steps are fused together to create a walk of M + N + 1 steps. We now 
define the observable of interest to be the indicator function defined as follows: 
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Figure 1: Concatenation of two walks on the square lattice. On the left the indicator function 
B{uii,ui2) = 1. while on the right B{uii,uj2) = 0. 



See Fig. 1 for two examples of concatenation. 

The more common definition for concatenation has the root points for the two walks placed at the 
origin. We use an alternate definition because it is straightforward to calculate the indicator function 
using our SAW-tree implementation [15]. 

If we let O be the coordination number of the lattice (0 = 6 for the simple cubic lattice), we then 
have 



Bm,n = Mean value of B{u}i,uj2) over all pairs of M and A'' step walks, 

= {B{'jJl,'^2))\ui\=M,\uj2\=N, 

B{uii,Ul2), 



— E 

^^^^ \u„\=M,\u„\=N 
Cm+N+1 



(3) 
(4) 

(5) 
(6) 



The longest walks which have been exactly enumerated on the simple cubic lattice have 36 steps [10], 
and we can recursively exploit this fact. For convenience we define 



Bjv = ^Bm N — 



N,N 



and so 



C2N+I 



Cur = Bl3Cr^ = BysBggCgg, 

2 54 8 
73-"36''36' 



C295 — -Bl47-B73-B3gcL, 



(7) 



(8) 
(9) 
(10) 



C38797311 = -Bl9398655-D9699327 " " ' -^36 C36 • 



(11) 



Thus, estimates for Bjy can be mapped to estimates of the number of walks cn- We can then use 
equation (1) to estimate /x: 



l/N 

UN = Cj!f 



log ^lN 



1 

N 
log/x 



log Cn 



(7-l)logiV^log^^^^^_^,_,^ 



N 



N 



(12) 
(13) 

(14) 



Corrections to scaling vanish with increasing N, and estimates for ji^ approach /i. 

Taking the logarithm of each side of equations (8)-(l 1), one can see that the contribution of the cae 
term remains approximately constant, but the addition of higher order terms successively eliminate the 
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higher order corrections. In particular, 



logM38797311 — 



1 



log -Bi9398655 + 



2 



log -B9699327 H 



38797311 



38797311 
220 



••• + 



38797311 



log-Bse + 



38797311 



l0gC36 



(15) 



The approach described here may be thought of as a "divide-and-conquer" algorithm, where a long 
SAW is successively spHt into halves. This is in stark contrast to typical growth algorithms such as 
PERM, where SAW (and other combinatorial objects) are incrementally built up step by step. 

2.3 Scale-free moves 

In order to accurately estimate ^ from equation (15), we must find an efficient way to estimate Bjy. We 
estimate by sampling pairs of SAW of length via the pivot algorithm, and then i?Ar is the time 
average of r2iJ(wi, W2). The observable B is not a global observable in the same sense as, for example, 
the square end-to-end distance: it clearly depends strongly on the details of the structure of each walk 
close to the concatenation joint. 

We now present a simple yet subtle argument to show that if we naively sample pivot sites uniformly 
at random, then Tint for B will be 0{N). We will assume throughout that we are considering pairs of 
walks of length N. 

First, let us define zero atmosphere SAW as those self-avoiding walks for which one of the ends has 
all neighboring sites occupied. It is well known that zero atmosphere walks have positive density in the 
set of all walks (see e.g. [16]). We denote a zero atmosphere SAW as "minimal" if, starting from the end, 
we visit all of the neighbors of the end in the fewest possible number of steps. Minimal zero atmosphere 
walks also have positive density in the set of SAW. E.g. for the square lattice, the density of minimal 
zero atmosphere SAW which start with the seven steps in Fig. 2 is bounded below as the SAW length 
iV -J> 00. 



Our ensemble is pairs of SAW, each of fixed length. Suppose we were to sample pivot sites uniformly 
at random, so generating a Markov chain. Assume we have equilibrated the Markov chain so that we 
are guaranteed to be sampling from the equilibrium distribution. If we were to choose a random time in 
the Markov chain, the probability of choosing a minimal zero atmosphere walk is then 0(1). However, 
the probability that the next pivot site chosen could change the value of the atmosphere is 0{1/N). 
Therefore, in this case B will, on average, remain zero for 0{N) time steps in the Markov chain. For the 
observable B, the contribution of zero atmosphere walks ensures that it must take time 0{N) to achieve 
an essentially new configuration. Thus, Tuii(B) = 0{N) when pivots are sampled uniformly at random. 
Note that this effect is actually quite subtle, as although zero atmosphere walks have positive density, in 
practice this density is small. Thus the contribution of these configurations to rjnt (^) is small in practice 
until N is of the order of thousands or tens of thousands. 

However, it is possible to dramatically improve the integrated autocorrelation time for B, and hence 
the accuracy of our estimate of B^. The key point is that the concatenation operation introduces a new, 
important length scale into the system, namely the distance from the concatenation joint to internal sites 
of the walk. B depends strongly on the structure of the walk according to this distance. We make the 
following conjecture. 




Figure 2: Minimal trapped walk of seven steps on the square lattice (solid Une) with a possible extension 
(dashed Une). 
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Conjecture 1 Suppose we have an observable for a polymer system that depends on a single internal 
distance, L. Then the integrated autocorrelation time for this observable is of the same order as the time 
it takes to make successful pivots at all length scales with respect to this distance. 

To be concrete, if L is the distance from an internal site to the concatenation joint, then we believe that an 
essentially new configuration with respect to B is obtained once pivots have been made at length scales 
Loforderl,2,4,8, 16, ••• ,7V. 

By choosing pivot sites uniformly at random with respect to log L, we therefore expect that there 
is only at most a log N penalty for the integrated autocorrelation time for as compared to a global 
observable, i.e. Tint (B) = 0{NP log N). N.B., since the CPU time per attempted pivot for the SAW-tree 
implementation is 0(log N) [15], this means that in CPU units n^tiB) = 0{NP log^ N). 



2.4 Experimental design 

To estimate p we must calculate each of the terms in equation (15). We do so by running separate Monte 
Carlo simulations for pairs of walks of length 36, 73, • • • , 19398655, in order to calculate B^. Since it 
takes CPU time 0(log N) to make a pivot attempt, and CPU time 0(log N) to calculate B, we choose 
to sample B for every time step in the Markov chain. The procedure we used was: 

1. Use the pseudo_dimerize procedure of [15] to generate two initial A'^-step SAW configurations. 

2. Initialize Markov chain by performing at least 20A'^ successful pivots on each SAW. Pivot sites are 
sampled uniformly at random. The stopping criterion must be based on the number of attempted 
pivots so as not to introduce bias. 

Our sampling procedure for B is then: 

1 . Select one of the two walks uniformly at random. 

2. Select a pivot site on this walk by generating a pseudorandom number x between and log N, and 

let pivot site j = [e^J . 

3. Attempt pivot move, update walk if result is self-avoiding. 

4. Randomly pivot each of the walks around their root points. These pivots are always successful. 

5. Calculate B{oJi ,oj2), and update our estimate of Bj^. 

6. Repeat. 

Our goal is to optimally partition CPU time amongst the terms in equation (15), in order to minimize 
the overall error in our estimate of p. The terms in equation (15) approach log-Bjv for large N. We 
have 

llogBAr^^logiV + 0(l/7V). (17) 

The 1/7V factor on the right hand side of the above equation dominates the increase in integrated auto- 
correlation time in CPU units for B. Therefore if we were to invest the same CPU time in each term of 
equation (15), the contributions to the error would dinninish with increasing A^! 

To minimize overall statistical error we now perform a short test run of CPU time for each length, 
determining the constants in 

(IS, 
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We show these measured values of a y in Fig. 3. However, we can also express a in terms of the variance 
of B and the integrated autocorrelation time of the algorithm. Assuming that Conjecture 1 is correct, 
modulo logarithmic factors we obtain the following expression for qn'- 

aN ~ 7V-l+(P+7-l)/2 ~ 7V-0-87. (19) 

In Fig. 3, it is clear that ajv decays as a power law with N, as expected. By inspection, ajv foUows 
the predicted power law behavior quite closely, and thus Fig. 3 provides strong numerical support for 
Conjecture 1. 




10^ 10^ 10^ 10^ 10^ 10^ 

N 



Figure 3: Measured values of ajv, which measures the expected error of contributions to equation (15), in 
units of \/ seconds. A hne of slope (— 1 + (p + 7— 1)/2) « —0.87 is included in the plot for comparison. 

We then fix the total running time for our computer experiment at t. The (statistical) square error in 
our estimate for jj, is then 

(7^ = ^ subject io t = ^^ti. (20) 

The optimal choice of ti to minimize is then 

ti = ^t, (21) 



and the optimal value for the error is 



(22) 



In practice, we did not rigorously apply this prescription to the longest walks, and instead spent at 
minimum 1 % of the CPU time at each length. 

Almost all of the computational effort is spent on the B^q and -B73 terms in equation (15). The 
higher order terms reduce the corrections to scaUng, and essentially eliminate the systematic error in our 
estimate of fx. 



3 Results and Analysis 

The analysis for this computer experiment is remarkably simple. It is an extremely rare example of a 
problem in lattice statistical mechanics for which we have strong evidence that systematic errors are 
negUgible. Hence the confidence intervals we report are purely statistical. 
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We ran the computer experiment for a total of 60 000 CPU hours on SunFire X4600M2 machines 
with 2.3GHz AMD Opteron CPUs. 

In Table 1 we report our estimates for Bn, and thence our estimates for cjv from equations (8)-(ll). 
Note that the estimates for are highly correlated. The error in the mantissa is given in the final column; 
for example, from Table 1 we estimate that C38797311 = 6.6 x io260i8276^ ^jjg confidence interval 
of the mantissa being (5.3, 8.2). This is a direct estimate from our values: it is not an extrapolation, 
and the reported error is purely statistical. As a technical aside, the error estimates for Bn in Table 1 are 
approximately constant for N > 151551 because we invested the same percentage of CPU time in each 
of these cases. 



N 



(Ar-l)/2 



Cat 



Cjv mantissa interval 



73 


2.47267030(65) 


2.139271x10^9 


(2, 


.139270, 


2.139271) 


147 


2.20753977(91) 


1.010276x10^9 


(1 


.010275, 


1.010277) 


295 


1.9740142(14) 


2.014793x10198 


(2. 


.014790, 


2.014796) 


591 


1.7668271(18) 


7.172241x10^96 


(7, 


.172218, 


7.172264) 


1183 


1.5823991(25) 


8.140025x10^93 


(8, 


.139971, 


8.140078) 


2367 


1.4178577(36) 


9.394724 xlQi^^^ 


(9. 


.394599, 


9.394850) 


4735 


1.2708081(58) 


1.121626x1031^'' 


(1. 


.121595, 


1.121656) 


9471 


1.1392521(81) 


1.433230x10^3^2 


(1. 


.433151, 


1.433308) 


18943 


1.0214669(91) 


2.098243x101=^^°* 


(2. 


.098013, 


2.098474) 


37887 


0.9159517(92) 


4.032592x10^5408 


(4, 


.031706, 


4.033477) 


75775 


0.8214372(97) 


1.335804x105°^!'' 


(1. 


.335217, 


1.336391) 


151551 


0.736643(10) 


1.314444xl0ioi'534 


(1. 


.313290, 


1.315600) 


303103 


0.660651(10) 


1.141449x10203268 


(1. 


.139445, 


1.143457) 


606207 


0.592531(11) 


7.720126 xlO-^oesss 


(7, 


.693038, 


7.747310) 


1212415 


0.531449(11) 


3.167451x10™! 


(3, 


.145262, 


3.189797) 


2424831 


0.476654(11) 


4.782146x101626142 


(4, 


.715379, 


4.849858) 


4849663 


0.427497(11) 


9.776394x103252284 


(9, 


.505309, 


1.005521) 


9699327 


0.383408(12) 


3.664531x106504569 


(3. 


.464124, 


3.876531) 


19398655 


0.343919(12) 


4.618409x1013009138 


(4, 


.127077, 


5.168235) 


38797311 


0.308455(11) 


6.579250x1026018276 


(5, 


.253839, 


8.239029) 



Table 1: Estimates of B^ and cjv with statistical errors. 

In our analysis for /x we utilize an estimate for the critical exponent 7 from a Monte Carlo computer 
experiment [17]: 7 = 1.15696(1). In addition, we utilize the estimate of the critical amplitude A = 
1.215(2) from [4]. We do this by setting 7* = 1.15696, A* = 1.215, and forming the improved 
estimates 

, (7*-l)logiV log A* 
log iiff = log UN - -j^ ^ • (23) 

We denote the errors in the utilized estimates as = 0.00001 and a a = 0.002. In the linait of large N, 
IJ,*j^ will then have the following contributions to the systematic error: 

The Ai term comes from the leading order correction in equation (14). From [14] we have Ai = 
0.528(12). The constant of this term is indeterminate, but we will see that it cannot be so large so as to 
interfere with our estimates. 

Our estimates for ji are collected in Table 2. The /x^ estimates rapidly converge with increasing N, 
which indicates that for the largest values of N systematic errors are negligible. We can also see from 
the table that the statistical error, a{ii*j^), is dominated by the low order terms. Finally, it is clear that the 
contributions from the errors of the 7* and A* terms are much smaller than the statistical error for large 
N. 
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One additional point is that for the largest values of N, N~'^'-~'^ is of the order of 10^^^. In principle, 
this term could have a large constant and result in a large and unknown systematic error. In practice, 
because of the smooth convergence of our estimates we know that the constant cannot be large, and 
hence contributions from this term to /x^ are negligible for large N. 



N 








\xa^\og N/N 


Hua/{AN) 


iV-^i" 


1 


73 


4.68373253707 


1.70x10" 


8 


2.79x10" 


■6 


1.07x10" 


-4 


1.60x 10" 


-3 


147 


4.68392658487 


2.13x 10~ 


■8 


1.60x10" 


-6 


5.28x10- 


-5 


5.61 X 10" 


-4 


295 


4.68400034315 


2.40 X 10~ 


■8 


9.06x10" 


-7 


2.62x10" 


-5 


1.97x10" 


-4 


591 


4.68402683289 


2.53x10" 


8 


5.07x10" 


■7 


1.31x10" 


-5 


6.96x10" 


-5 


1183 


A a A r\'~\ c r\ A '—1 

4.68403589477 


2.60x 10~ 


8 


2.80x10" 


-7 


6.52x10" 


-6 


2.46x10" 


-5 


2367 


4.68403883775 


2.64x 10~ 


■8 


1.54x10" 


-7 


3.26x10" 


-6 


8.68x10 


-6 


4735 


4.68403971655 


2.68 X 10 


8 


8.37x10" 


-8 


1.63x10" 


-6 


3.07x10" 


-6 






Z. /W X iU 


8 


4.53x10" 


-8 


8.14x10- 


-7 


1 .Uo X lU 


-6 


18943 


4.68403997588 


2.71x10" 


■8 


2.44x10" 


-8 


4.07x10" 


-7 


3.84x10" 


-7 


37887 


4.68403996593 


2.71x10" 


8 


1.30x10" 


■8 


2.04x10" 


-7 


1.36x10" 


-7 


75775 


4.68403995443 


2.71x10- 


8 


6.95x10- 


9 


1.02x10- 


-7 


4.79x10" 


-8 


151551 


4.68403994395 


2.71x10" 


■8 


3.69x10" 


-9 


5.09x10" 


-8 


1.69x10" 


-8 


303103 


4.68403993749 


2.72x10" 


8 


1.95x10" 


9 


2.54x10" 


-8 


5.99x10" 


-9 


606207 


4.68403993406 


2.72x10- 


8 


1.03x10- 


9 


1.27x10- 


-8 


2.12x10" 


-9 


1212415 


4.68403993235 


2.72x10" 


■8 


5.41x10" 


-10 


6.36x10" 


-9 


7.49x10" 


-10 


2424831 


4.68403993145 


2.72x10" 


8 


2.84x10" 


■10 


3.18x10" 


9 


2.65x10" 


-10 


4849663 


4.68403993096 


2.72x10" 


8 


1.49x10" 


-10 


1.59x10- 


9 


9.36x10" 


-11 


9699327 


4.68403993069 


2.72x10" 


■8 


7.77x10" 


-11 


7.95x10" 


-10 


3.31x10" 


-11 


19398655 


4.68403993058 


2.72x10" 


8 


4.05x10" 


-11 


3.97x10" 


-10 


1.17x10" 


-11 


38797311 


4.68403993052 


2.72x10" 


8 


2.11x10" 


-11 


1.99x10" 


-10 


4.14x10" 


-12 



Table 2: Estimates of with statistical error fT(/ti^), and contributions to the systematic error. 

We thus conclude that the estimate /i387973ii has negligible systematic error, and hence adopt this as 
our best estimate for /i. Our final estimate is /i = 4.684 039 931(27). 

Note, we could have avoided the use of previous estimates of 7 and A, had the calculation of Bj^ 
been extended to larger N. This was not done because for N of the order of 100 million or so, both 
memory management and initialization time become significant but not insurmountable issues for the 
simulation of SAW using the SAW-tree implementation [15]. 

4 Discussion 

As noted in the introduction, for the calculation of ji the approach which is most competitive with the 
algorithm presented in this paper is PERM [9], where the estimate n = 4.684038 6(11) was obtained. 
Our error bar is approximately 40 times smaller, which is clearly a significant improvement upon the 
previous state of the art. Other approaches to the calculation of ji worth noting are the method of atmo- 
spheres [18], and the Berretti-Sokal algorithm [19]. 

We note in passing that the method of atmospheres could be combined with the pivot algorithm and 
scale-free moves to obtain an accurate estimate for ji. We will not go into any depth, but the method of 
atmospheres corresponds to estimating 

CN+K _ 

cnCk Afi^CK Ck 

for small, fixed K, and in the limit N ^ 00. From this expression one can then estimate fi once 
corrections-to-scaling have been taken into account. Despite being more accurate than previous methods, 
it is an order of magnitude less accurate than the method described here. This is because the mean CPU 
time per pivot attempt is 0{logN) for the SAW-tree implementation. For the atmospheric sampling 



8 



method, the dominant error comes from sampling walks in the large A'^ limit, while for the method 
described in this paper the dominant error term comes from samphng short walks (in our case, with 
N = 36). 

On the topic of approximation enumeration of SAW beyond the limit of exact enumeration, there have 
been a number of papers in recent years. Approaches include incomplete enumeration [20], flatPERM 
and flatGARM [21], stochastic enumeration [22], and the multicanonical Monte Carlo method [23]. 
The relative advantage of our approach is significant for small N, e.g Shirai and Kikuchi [23] obtained 
C256 = 6.2(4) X 10^°^ for the square lattice, while for comparison we found C295 = 2.014793(3) x 10^^* 
on the simple cubic lattice. For larger N, the relative advantage of our method increases, since to generate 
a SAW using an incremental growth method takes CPU time at least 0{N). This factor of N becomes 
prohibitively large when N is of the order of millions. 

It is not clear to us if our approach could be adapted to other approximate enumeration problems, 
or to estimations of the free energy for other models in statistical mechanics. The general principles of 
"divide-and-conquer" and the use of global moves in the canonical ensemble may be of wider use, or it 
may be that SAW is a particularly favorable model. 

We consider Fig. 3 to be strong evidence in favor of the correctness of Conjecture 1. We therefore 
expect that the use of scale- free moves for the simulation of polymers will prove useful in other contexts 
where there are additional length scales. For example, in the cases of star polymers or confined polymers. 
We will explore this idea further in a future paper where we wiU also derive an estimate of the critical 
exponent 7 [17]. 

In future, our implementation of the SAW-tree could be optimized for the non-uniform selection 
of pivot sites according to our scale-free prescription. In particular, there is no reason a pivot being 
performed near the end of a walk should take mean CPU time 0(log N). It is possible to arrange the 
binary tree data structure so that this operation would take time 0(1). One natural way of doing this 
would be to use a splay tree [24], which would dynamically adjust to form an optimal tree structure for 
any choice of pivot site sampling distribution. 

We could also obtain a constant factor improvement, if it were possible to efficiently forbid configu- 
rations with immediate returns at the concatenation joint. 

Finally, it is certainly possible to apply this approach to other lattices. Unfortunately, in the case of 
the square lattice the finite lattice method enumerations of polygons provide estimates for /x [8] which are 
approximately 2 orders of magnitude more accurate than our method. However, for three-dimensional 
lattices such as the body centered cubic lattice and the face centered cubic lattice, our method will allow 
for much more accurate calculations of n than are currently available. 

5 Conclusion 

We have apphed the pivot algorithm to calculate the coimective constant for self-avoiding walks on 
the simple cubic lattice, obtaining /i = 4.684 039 931(27). Our approach may also be used to derive 
extremely accurate estimates for the number of self-avoiding walks. The power of our approach derives 
from the apphcation of an efficient global move (the pivot algorithm), use of an observable which is 
calculated through a divide-and-conquer approach, and from the apphcation of scale-free moves. We 
hope that these key ideas may prove useful in other contexts. 
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